Fractional Calculus and the Evolution of Fractal Phenomena 
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It is argued that the evolution of complex phenomena ought to be described by fractional, 
differential, stochastic equations whose solutions have scaling properties and are therefore random, 
fractal functions. To support this argument we demonstrate that the fractional derivative (integral) 
of a generalized Weierstrass function (GWF) is another fractal function with a greater (lesser) fractal 
dimension. We also determine that the GWF is a solution to such a fractional differential stochastic 
equation of motion. 
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I. INTRODUCTION 

Significant changes in our knowledge of how to analyze 
nonlinear static and dynamical phenomena in the phys- 
ical and biological sciences have occurred over the past 
thirty years |l]-[|]. In the physical sciences methodology 
has moved away from complete reliance on the tools of 
linear, analytic, quantitative mathematical physics to- 
wards a combination of nonlinear, numerical and quali- 
tative techniques ||. Not only have many of the linear 
analytical approaches often proved to be inadequate, but 
the entrenched geometry of Euclid, classically used to de- 
scribe natural phenomena, has not always been adequate 
to the task. In the 1960s Mandelbrot began discussing a 
new geometry of nature, one that embraces the irregular 
shapes of objects such as coastlines, lighting bolts, clouds 
and molecular trajectories. A common feature of these 
objects, which Mandelbrot called fractals, is that their 
boundaries are so irregular that it is not easy to under- 
stand how to apply simple metrical ideas and operations 
to them Towards this end we shall consider some 

of the metric peculiarities of a few unusual mathematical 
objects, fractal functions, and discuss the possible physi- 
cal implications of their evolution in time in terms of the 
phenomena they are used to model. 

One of the defining properties of a fractal function 
is that it does not possess a characteristic scale length 
and consequently its derivatives diverge. The Weierstrass 
function was the first example of a function that is contin- 
uous everywhere, but is nowhere differentiable. In 1926 
Richardson || suggested that this function, or one shar- 
ing its non-analytic properties, be used to describe the 
velocity field of the atmosphere, because of the appar- 
ent impossibility of modeling the irregular flow of the 
wind using differential equations. It was specifically this 
property of non-differentiability that Richardson believed 
captured the essential features of turbulence. This view 
of turbulence, rather than becoming quaint and out of 
date, has been demonstrated to be quite modern and 
in fact anticipated the introduction of fractals into the 
description of complex phenomena . 



During the same period Levy ]7j was working to es- 
tablish the most general properties necessary for a pro- 
cess to violate the then accepted form of the Central 
Limit Theorem and still converge to a limit distribu- 
tion. He was quite successful, establishing the class of 
infinitely divisible distributions, which as its name im- 
plies concerns processes whose statistical properties per- 
sist at each level of aggregation of the data, and are to- 
day called a-stable Levy processes. A deep connection 
between Levy stable processes and the Weierstrass func- 
tion was established using random walk concepts |^ [n]] 
and subsequently used to understand turbulent fluid flow 

§■ 

It is not only macroscopic phenomena such as turbu- 
lence that suffer from the problem of not being describ- 
able by analytic functions, however. The theory of Brow- 
nian motion as formulated in 1908 by Langevin Jll]] has 
the form 



du (t) 
dt 



-Xu (t) + F (t) 



(1) 



where u (t) is the velocity of the Brownian particle, A 
is the dissipation rate and F (t) is the random influence 
of the bath of lighter particles (molecular impacts) on 
the Brownian particle. In 1909 Perrin |l^| observed that 
the path of a Brownian particle, as seen through his mi- 
croscope, is continuous but not differentiable. Therefore 
such a path can not be described by an analytic function. 
In 1942 Doob []l3| proved that the velocity of the Brown- 
ian particle is discontinuous for the Ornstein-Uhlenbeck 
process Jl^] described by ([!]) and discusses how one con- 
structs a stochastic differential equation to replace the 
more familiar differential equations of physics up to that 
time. In particular he replaced ([!]) with 



du (t) = -Ait (t) dt + dB (t) 



(2) 



where dB (t) is the increment of a Wiener process. This 
new mathematics has become so accepted that we now 
write the Langevin equation (^), but interpret it as the 
Doob equation (|2|). A process appears random on the 
macroscopic level due to the separation of time scales in 
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the microscopic and macroscopic processes. When such a 
separation exists the Langevin equation, using the inter- 
pretation of Doob, adequately describes the dynamics of 
the physical phenomenon. On the other hand, when this 
separation of time scales does not exist, ordinary statis- 
tical physics is no longer adequate to describe the phe- 
nomenon, as discussed by Grigolini et al. Jl5[ |. In partic- 
ular, a lack of time-scale separation induces a fractional, 
stochastic, differential equation on the macroscopic level 

§- 

What the Weierstrass and other non-analytic functions 
and the Levy distribution have in common is the lack of 
a differential description of the time evolution of the phe- 
nomena which they are intended to model. It was only 
a decade ago that it was recognized by Shlesinger et al. 
H that Levy flights can be used to describe the velocity 
field of turbulent flow. The key to this understanding 
is that the equation of evolution of the probability den- 
sity is a fractional diffusion equation whose solution is a 
truncated Levy distribution fl6fl . 

Fractional diffusion equations have been used to model 
the evolution of stochastic phenomena with long-time 
memory, that is, phenomena with correlations that decay 
as inverse power laws rather than exponentially in time 
16 18 1 . It has been known for quite some time in the 



economics literature that such statistical properties can 
be successfully described using finite difference equations 
in which the finite difference is of fractional order, see 
Hosking Jll],|2(| . The continuum limit of such fractional 
difference stochastic equations are fractional differential 
stochastic equations, 



d^u jt) 



F(t) 



(3) 



where here the mathematics has been generalized to 
the fractional calculus to accommodate the non-analytic 
properties of the underlying phenomena. The evolution 
of the probability density associated with the velocity 
in (^) is a fractional partial differential equation in the 
phase space for the phenomenon, see for example West 
and Grigolini [^lf . 

Herein we suggest that the dynamics of complex phe- 
nomena, described by fractal functions, can be expressed 
in terms of fractional differential equations of motion. 
One approach to doing this has been adopted by Non- 
nenmacher and colleagues B2 ], who generalize the tradi- 
tional models of viscoelastic phenomena using fractional 
initial value problems. The solutions to their fractional 
equations of motion yield physical observables, such as 
the stress relaxation, that is in excellent agreement with 
experiment. Others have generalized the analysis of wave 
phenomena, starting with the diffractals of Berry [ p3[ , up 
to and including the fractional wave equation of Schnei- 
der and Wyss Q. 

Our goal here is more modest in that we do not attempt 
to describe any particular phenomenon in complete de- 



tail. On the other hand it is more ambitious in that we 
demonstrate the generic result that fractal functions have 
fractional derivatives and therefore complex phenomena 
having a fractal dimension are more reasonably modeled 
using fractional equations of motion than they are using 
ordinary differential equations of motion such as given by 
Newton's laws. 

Even sharing the universally accepted view according 
to which Newton's laws apply to the fundamental con- 
stituents of complex phenomena |p5| , this microscopic 
level of description is often not the natural one for the 
phenomena being studied. For example, one rarely uses 
particle dynamics to describe the viscous fluid flow in tur- 
bulence, instead one more often uses the Navier-Stokes 
equation to describe the behavior of the velocity field 
pfjj. The familiar, ordinary differential equations of mo- 
tion for individual particles need to be replaced with frac- 
tional differential and integral equations for the appropri- 
ate field variables. We show, using a generalized Weier- 
strass function (GWF), that if said function has a fractal 
dimension D then the /3-fractional derivative yields a new 
function of fractal dimension D + (3. In a similar way the 
/3-fractional integral of the generalized Weierstrass func- 
tion yields a new function of fractal dimension D — (3. 
Using these properties we can write the fractional dif- 
ferential equations of motion for the complex process of 
interest. 



II. GENERALIZED WEIERSTRASS FUNCTION 

The Weierstrass function was the first exemplar of a 
function that is continuous and non-differentiable. A 
number of generalizations of this fractal function have 
been discussed in the literature Jj]|2],|2"7|] , the most gen- 
eral including random behavior. We refer to the latter as 
the generalized Weierstrass function (GWF) defined by 
Berry and Lewis [p| as the series: 



W(t) 



E 

n— — oo 



(i 



y(2-D)n 



(4) 



where 1 < D < 2, 7 > 1 and <f> n is an arbitrary phase. 
Berry and Lewis show using a combination of numerics 
and analysis that for these values of the parameters the 
GWF is continuous but not differentiable, has no charac- 
teristic scale, and almost certainly has a fractal dimen- 
sion D [^9f . This argument was given mathematical rigor 
by Mauldin and Williams |30) somewhat later. Further, 
if the set of phases {4> n } are uniformly distributed on the 
interval (0, 2ir) then W(t) is a stochastic fractal function. 
The properties of a modified version of (Q) was considered 
by Falconer p7| . 

It is well known that (Q) has a fractal graph, that is, the 
trail left by the GWF in (X, t)-space is a fractal curve. 
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The fractal nature of such graphs are manifest in the 
power-law behavior of the correlation between measure- 
ments separated by a time interval r p^] . Before we con- 
sider the increments in the GWF, that is, the difference 
in measurements separated by a time r, let us examine 
the scaling properties of (0) 



W(jt) = -f 2 - D W(t) 



(5) 



which is obtained by relabeling the series index. The 
most general solution to the scaling equation (0) is pLpf 



W (t) = A (t) t c 



(6) 



where A (t) is a function periodic in the logarithm of t, 
with period hry, and the power-law index is given by 



a = 2-D. 



(7) 



Thus, we see that the dominant behavior of the GWF is 
a power-law growth in time, and the power-law index is 
determined by the fractional dimension, D. 

It is possible to connect the expression for the frac- 
tional dimension, (ffl), with the slope of the correlation 
function of the increments in the GWF. Let us define 
such a correlation function in the following way 



C(r) = (\AW(t,r)\ 2 )t 



(8) 



where the increments of the GWF are defined by 

AW(t,r) = W{t + r)-W{t) (9) 

and the average, indicated by the brackets with a <fr sub- 
script in (|g) is taken over an ensemble of realizations 
of the phases {(/>} uniformly distributed on the interval 
(0, 2-71"). We refer to the function (||) as the incremental 
GWF or IGWF. Inserting the IGWF 



AW(*,t) = 



71= — OO 



7 



-(2-D)n 



(1 - e !7 T )e 



(10) 



into the phase average in (^|) and carrying out the phase 
average 



«■■■».= II jf £<■■■> 



("I 



TI— — OQ 



which for a uniform distribution of phases yields 

(e*<*»-*»'>) =S n , n >, (12) 
the correlation function becomes 

OO 

C(t)=2 7~ 2(2 ~ D)ii (1-cos 7 "t). (13) 

Tl— — OO 

Note, since the correlation of the IGWF are independent 
of the initial time t, these increments are a realization of 
a stationary stochastic process. 



The dominant behavior of the IGWF correlation func- 
tion ( |l3j ) is determined by the solution to the exact scal- 
ing relation 



C( 7 t)= 7 2(2 ~ D) C(t) 



(14) 



obtained from (|T^) by again relabeling the series. The 
most general solution to ( 14 ) is 



C(t) = A(t)t 2 



(15) 



where as before A(t) is a function periodic in the loga- 
rithm of its argument and the power-law index is given 
by ([?]) [0||. Note that a plays the same role here as 
the Hurst exponent, H, does in random walk processes. 
In fact if a = H = 1/2, then the correlation function 
increases linearly with time, so that the IGWF, AW, 
would be a normal diffusive process with Gaussian statis- 
tics that is stationary in time. 



III. FRACTIONAL CALCULUS AND THE IGWF 

We are interested in both the fractional integral and 
the fractional derivative of the generalized Weierstrass 
function. However, in applying the fractional calculus, 
we take note of the fact that the GWF is not a station- 
ary stochastic function. It is the IGWF that is stationary. 
Therefore we apply the fractional calculus to the station- 
ary IGWF and thereby avoid some technical difficulties. 



A. Fractional integral of the IGWF 

Let us introduce the Riemann-Liouville definition of a 
fractional integral of order (3 of the GWF: 



4 W(t')dt' 



r(/3) V-oo (t-t'y-P" 



(16) 



where < (3 < 1. We have for the fractional integral of 
the IGWF: 

^^- m Lw*w** wv ' T) - (17) 

so that inserting ( |l0| ) into (jDj) yields 



1 ^ e* 



T{P) f-' 7 (2-D)n 



(1 - e 17 T ) 



71 — — OO 

* dt' 



The integral between the curly braces yields 



e 17 "*' S> (18) 



3 



T(J3)^-e 

and therefore 
AW(-V(t,T) = e-™V 2 



t/3/2 



(19) 



(2-D+0)n 



n— — oo 



7 



(l-e l7 " T ). 



(20) 



Thus, we see that the fractional integral of the IGWF 
has the same form as the original IGWF. The difference 



between (20) and (|10|), up to an overall phase, is that 
D — > D — 3. What is the proper interpretation of this 
shifting of the parameter value? 

To answer this question, we need to go back to the 
GWF. We shall address this issue in the next section. 
For the time being, we limit ourselves to noticing that 
the correlation function related to ( |2C| ) is given by 

oo 

C(t) = 2 7~ 2(2 ~ D+/3) "(1-cos7™t) (21) 

n=— oo 



and therefore, if fl20|) corresponds to some difference of 
properly defined GWFs, by virtue of the scaling index 
(0), these new GWFs would have a new fractional di- 
mension given by D' = D — fi. In the next section we 
shall make this argument rigorous. 



B. Fractional derivative of the IGWF 

The calculations in the case of the fractional deriva- 
tive of the IGWF are similar to those carried out in the 
case of the fractional integral of the IGWF. Let us con- 
sider the Ricmann-Liouvillc definition of the 8 fractional 
derivative of the GWF: 



V^W(t) EE 



1 



* W{t')dt' 

r(i - 0) dt (t - t' f 



(22) 



where again < j3 < 1. The fractional derivative of the 
IGWF result in 

(23) 

The expression (|18|) is now replaced with 



1 oo 

A^(*,r) = -4 £ 



J4n 



TIB) ^ 7 (2-£>)n 

v ! n — — OO ' 



dt' 



dt ./_«, (t - ty-e 



(1 - e* 7 T ) 



Jj n t' 



(24) 



where the time derivative of the integral in curly braces 
is 



fT(l -3) 7 /3« e i7 n t e -^(^-l)/2 

and therefore 



n— — oo 



e »(</>„+7"t) 
7 (2-D-/3)n 



(1 



,17 r\ 



(25) 



(26) 



Thus, we see that the fractional derivative of the IGWF 
has the same form as the original IGWF. The difference 
between (|2^) and (|l0|), up to an overall phase, is that 
D -> D + /3. 

Again after observing that the correlation function re- 
lated to ( p6| ) is 

oo 

C*(t) = 2 ^ 7 -2(2-D-/J)n ^ _ COS7 « r ) (27) 
n— — oo 

we make the plausible conjecture that D + 8 corresponds 
to the new fractional dimension, D", of some properly 
defined GWF. Again this is made rigorous in the next 
section. 



IV. WHAT ABOUT THE GWF? 

We now want to use the fractional integral and frac- 
tional derivative of the IGWF to determine these same 
operations of the GWF. To accomplish this we assume 
that the IGWF is a given, as are its fractional integral 
and fractional derivative, but the GWF remains to be 
determined. Using the RHS of (|) and @ we obtain 



W (t + t) - W (t) = 



71— — OO 



7 



(2-D)n 



E 



s i( 7 "(t+r)+0„) 
j(2-D)n 



from which we can define the function 

J,(i n t+<f> n ) 



f(t) = w(t)+ J2 



7 



(2-D)n 



(28) 



(29) 



where W (£) is assumed to be unknown, and ( |28| ) is re- 
placed with the condition 



fit- 



fit) 



(30) 



The series in (|29|) is divergent, so that in order to regu- 
larize the function we write 



iV 



f N (t) = w N (t)+ 



„i(7"t+<M 



-N 



7 



(2-D)n 



(31) 



where 



4 



/ (i) = lim f N (t) 

iv— >oo 

W(t) = lim W N (t). 
The constraint (|30| ) is now replaced with 



lim f N (t + r) = lim /jv (t) 

iV— >oo N^oo 



(32) 



(33) 



indicating that the regularized function is either periodic 
with period t or is a constant in the limit. 

We now wish to establish that knowing the IGWF 
uniquely determines the function W (i) , and this is the 
GWF. To accomplish this we use the constraint given by 
(|33|). We present the analysis for the case Jn {t) = An, 
where An is a constant, so that (J3l|) can be written 



W N (t) = A N 



N 



n=-N 



7 



(2-D)n 



(34) 



In order to determine the constant in ( |34| ) we impose 
an additional constraint on the function. Since we want 
the function W (t) to be a fractal we require that Wn (t) 
satisfy the scaling condition, 



lim W N ht) = lim bW N (t) 

N— >oo N^oo 



(35) 



where 6 is a constant. Note that we exclude the possibil- 
ity that the function Jn (t) is periodic, since this will not 
satisfy this additional requirement that the function also 
scales. Imposing this scaling constraint on (EmI) we have 



Wjv (7*) = A N - 7 



(2-D) 



N-l 

E 



e i( 7 "*+0„) 

j(2-D)n 



(36) 



so if we choose 6 = 7 
TV — > 00 limit if we also choose 



n=— W+l 

1 ° we can satisfy (|35|) in the 



.4 



N 



N 

E 

i=-N 



7 (2-D)n • 



where {0 n } is a set of arbitrary phases, because 



lim An = lim 7 



2-D 



A 



N- 



(37) 



(38) 



Thus, the divergences that required the regularization 
exactly cancel in ([56]) with the choice of constraint j37| ) 
and we obtain the function 



W(t) 



lim W N (t) 



1 



71 — — OO 



j(2-D)n 



(39) 



which is the GWF, where we have associated the set of 
phases with those used earlier in the GWF. 

This argument can also be applied to the fractional in- 
tegral of the IGWF given by (|2^) so that we also have 
for the fractional integral of the GWF 



7 



(2-D+/3)n 



(1-e^"*). (40) 



In the same way the fractional derivative of the IGWF 
given by ( p6|) implies that the fractional derivative of the 
GWF is 



W {m (t,T)=e' 



t/3/2 



E 



y(2-D-P)n 



(1-e 47 "*). (41) 



Thus, our earlier remarks regarding the fractional calcu- 
lus applied to the IGWF apply equally well to the GWF. 



V. CONCLUSIONS 

We have determined the fractional dimension of both 
fractional integral and derivative of the GWF. Our result 
reads in these two cases: 



Dim[I? ( - /3) IF] = Dim[VF] - /?, 
Dim[V w W] = Dim[PF] + (5. 



(42) 
(43) 



This result can be easily interpreted noticing that the 
fractional dimension gives information about the degree 
of irregularity of the function under analysis. We have 
demonstrated that carrying out a fractional integral of 
the GWF means decreasing its fractional dimension and 
therefore smoothing the process, while carrying out a 
fractional derivative means increasing the fractional di- 
mension and therefore making the process and its incre- 
ments more irregular. 

A related result was obtained by Kolwankar and Gan- 
gal I plj , ^ , but required the introduction of a local frac- 
tional derivative (LFD), that is, a fractional derivative 
defined such that its non-local character is removed. 
They find that the LFD of a Weierstrass function, that 
is, the imaginary part of (Q) with all cf> n = 0, exists up 
to "critical order" 2 — D and not so for orders between 
2 — D and 1, where D (1 < D < 2) is the box counting 
dimension of the graph of the function. It is possible to 
show that our result is in complete agreement with that 
of Kolwankar and Gangal 
consider the inequality 



To this aim, let us 



D T < D <D t 



(44) 



where Dt is the topological dimension and De is the em- 
bedding dimension JIJ]. The condition (Q) needs to be 
fulfilled by any function in order to be a fractal. In the 
specific case of the graph of a fractal function, Dt = 1 
and De = 2. The same condition (H) must also hold 
true for the fractional derivative (fractional integral) of 
the function, in this case the GWF, in order to preserve 
its fractal properties. Therefore, for the fractional deriva- 
tive of GWF, we have: 
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D"<2=>D + (3<2^I3<2-D. 



(45) 



This means that the generalized Weierstrass function is 
fractionally differentiable for all orders less than 2 — D, 
which is the result of Kolwankar and Gangal |3l],|32| . 

Moreover, the same kind of reasoning can be applied 
to the fractional integral of the GWF. In this case, the 
meaningful condition reads: 



D' >1=>D-(3>1^>(3<D-1. 



(46) 



This means that the generalized Weierstrass function is 
fractionally integrable for all orders less than D — 1. 

Finally, we go back to the main issue, that being, a 
possible equation of evolution for a complex system, ex- 
hibiting fractal behavior and perhaps representable by 
the increments of a generalized Weierstrass function. 

Consider the function 



/(*) 



A n (l 



J4>n 



(47) 



where {w n } is a set of frequencies, {4>n} is a set of ran- 
dom phase confined to the interval (0, 2tt), and A n is a 
real amplitude. Thus, / (t) is a stochastic function of 
time with the definite initial condition / (0) = 0. If we 
now use ( E^ ) as a driving force for a fractional stochastic 
equation we can write 



V {a) F{t) = / (t) 

as our equation of motion. The solution to 
mally given by the inverse equation 

F(t) = V { - a) f{t) 



(48) 
is for- 



(49) 



so that using the definition of the fractional integral, (16) 
we obtain the explicit form of the solution 



oo . 

*"(*)= E £(i 



J<t>'-n 



(50) 



where we have absorbed the overall phase into the new 
random phases {4>' n }- Therefore, if we choose the param- 
eter values u> n — 7™, A n — 1 and a = 2 — D, the solution 
to the fractional stochastic differential equation, (^), is 
the GWF. 

In addition if we choose the coefficients in the series 
representation of the random driving force by A n oc 1 /u> n 
then ( ^7| ) is a realization of a complex Brownian motion 
process, which is to say that the statistics are two dimen- 
sional Gaussian and the spectrum is an inverse square 
frequency. In this case if we also choose for the parame- 
ters in the solution to the fractional stochastic differential 
equation u) n = 7™ and a = 1 — D, then the solution is 
again a GWF. In addition, as pointed out by Mandel- 
brot |jj (page 390), the GWF with random coefficients is 
a good approximation to a fractional Brownian function. 
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